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I.   INTRODUCTION 

In  the  last  few  years  there  has  been  growing  evidence 
that  barotropic  instability  may  be  a  common  occurrence  at 
the  200-mb  level  in  the  tropics.   Krishnamurti  [1970,  1971] 
observed  a  peak  at  wave  numbers  6-8  in  the  wind  spectrum  at 
this  level  which  he  proposed  was  caused  by  barotropic  insta- 
bility.  A  subsequent  study  by  Kanamitsu,  et  al.  [1972]  showed 
that  for  the  region  between  15S  and  15N  waves  of  this  scale 
do  receive  energy  through  barotropic  interaction  with  the 
zonal  mean  flow  and  zonal  wave  number  1. 

Kuo  [1949]  examined  the  barotropic  stability  character- 
istics of  a  parallel  zonal  current  and  found  that  instability 
is  possible  only  if  the  gradient  of  the  absolute  vorticity 
changes  sign  somewhere  in  the  region.   The  200-mb  level  in 
the  tropics  is  the  outflow  level  of  the  giant  monsoon  circu- 
lation which  is  induced  by  differential  heating  (both  sensible 
and  latent)  due  to  land-sea  distribution.   The  basic  flow 
pattern  thus  varies  in  the  planetary  scale.   Detailed  obser- 
vations have  found  that  synoptic-scale  moving  disturbances 
at  this  level  tend  to  form  in  the  easterly  jet  south  of  the 
Tibetan  high  [Krishnamurti,  1970]  during  the  northern  hemis- 
phere summer.   This  is  a  region  of  large  vorticity  gradient 
where  the  necessary  condition  for  instability  is  occasionally 
satisfied  locally.   If  these  waves  were  the  result  of  baro- 
tropic instability,  they  would  extract  energy  from  the  mean 


flow  and  the  long  waves,  since  the  long  waves  and  the  mean 
flow  combine  to  give  large  vorticity  gradients  in  these 
regions.   The  energy  exchanges  computed  by  Kanamitsu  et  al . 
[1972]  are  consistent  with  this  view. 

In  this  thesis  the  barotropic  instability  of  an  upper 
level  easterly  jet  was  examined.   In  addition  to  the  calcu- 
lation of  the  growth  rates,  the  downward  propagation  of  energy 
was  also  investigated  to  examine  the  possibility  that  lower 
troposphere  tropical  disturbances  may  originate  at  the  upper 
levels.   Simple  CISK  type  heating  was  introduced  to  ascertain 
how  it  affects  the  growth  of  the  waves. 

This  study  uses  a  quas i-geostrophic  two-level  model.   The 
equations  are  linearized  and  solved  by  the  initial  value 
technique.   Growth  rate  and  wave  structure  are  extracted  after 
an  exponential  growth  is  established. 


II .   THE  FORECAST  EQUATIONS 

The  quasi-geostrophic  prediction  equations  were  used  as 
the  forecast  equations.   It  is  recognized  that  the  quasi- 
geostrophic  approximation  is  not  very  accurate  near  the 
Equator,  but  Matsuno  [1966]  has  shown  with  a  barotropic  model 
that  quasi-geostrophic  motions  are  identifiable  near  the 
Equator  if  the  wavelength  is  not  too  long.   Also,  in  a  study 
with  a  baroclinic  model,  Yamasaki  [1969]  has  found  that  the 
error  is  small  at  twenty  degrees  latitude. 

The  quasi-geostrophic  equations  are  simplified  through 
the  use  of  a  two-level  model.   The  basic  model  is  very  similar 
to  the  models  used  by  Williams  et  al.  [1971]  and  Robertson 
[1973].   The  equations  were  linearized  and  solved  numerically 
by  the  initial  value  technique.   The  growth  rates  and  wave 
structures  were  then  computed. 

The  simple  two- level  model  was  constructed  by  dividing 
the  atmosphere  into  four  layers  with  a  constant  pressure  dif- 
ferential of  Ap/2  (Fig.  1),  numbered  0  to  4  from  top  to  bottom 
The  boundary  condition  at  p  =  p   was  w   =0,  which  prohibits 
vertical  propagation  of  energy.   Charney  [1969]  has  shown  that 
such  energy  propagation  is  unlikely.   At  p  =  p   the  boundary 

condition  was  to  =  to   where  to   was  the  frictionally  induced  to 

e        e  J 

at  the  top  of  the  Ekman  layer.   Charney  and  Eliassen  [1949] 
derived  the  following  expression  for  w  . 
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WP  =  "P  g^Tr-l  /2    sin  2  a  c  (2.1) 

where  A   is  the  kinematic  eddy  viscosity,  a  is  the  surface 
inflow  angle,  and  j;  4  is  the  surface  geostrophic  vorticity. 

p  Level 

pt %  =  0 0 

ys    1 

Pt  +  AP a)2  2 

1  3 

po w  =  ^ 4 

Fig.  1.   Two-level  prediction  model 

It  has  been  shown  by  Holton  et  al .  [1971]  and  Chang  [1973] 
that  this  expression  is  reasonably  valid  as  far  equatorward 
as  ten  degrees,  but  may  be  invalid  close  to  the  Equator  south 
of  the  ITCZ.   The  surface  geostrophic  vorticity  was  approximated 
by 

U  -  £3  =  V2  ^  (2.2) 

The  quasi-geostrophic  vorticity  equation  in  pressure  coordin- 
ates is 


I^V2^  +   IKXVT  •  V(V2Y)  +  30-||   "  £0-||  =  0  (2.3) 

where  f   =  2  ft  sin  \b    .    6   =  2  Q  cos  ]U   /a,  w  =  dp/dt,  y  =  stream 
o  To    o  r0 

function,  ip   is  the  reference  latitude,  and  a  is  the  radius 
of  the  earth.   Observe  that  when  co  =  0,  (2.3)  becomes  the 


11 


barotropic  vorticity  equation  which  is  sufficient  to  describe 
the  barotropic  instability. 

Applying  this  equation  at  levels  1  and  3  gives 


&  V2^  +  fCXVf,  •  V(V^J  +  3Q  ^   -  fQ  ^  =  0,       (2.4) 


JL  V2^      +     pCXVY,    •    V(V2^ )    +   3      f^i-     -   f      ("*  -  "O    =    o. 
dt  3  r  3  v         3^  o     dx  OAP 


(2.5) 


If  it  is  desired  to  stop  at  some  intermediate  point,  say  the 
tropopause,  rather  than  to  include  the  entire  atmosphere, 
smaller  layers  could  be  used. 

Next,  consider  the  quasi-geostrophic  first  law  of  thermo- 
dynamics in  the  form 

^|I+  M„  .  V(|I)  ♦■£»  =  0  (2.6) 

where 

o  =  **1sl  [dTa  +  JL]  (2.6a) 

p^g  L  dz     C  J  '  v 

and  Ts  is  the  temperature  obtained  from  the  standard  atmosphere. 

The  heating  follows  the  formulation  of  Ooyama  [1964]  , 
Charney  and  Eliassen  [1964],  and  Kuo  [1965],  i.e.,  the  con- 
densation heating  is  proportional  to  the  convergence  of  moisture 
in  the  Ekman  layer.   The  form  used  was  given  by  Ogura  [1964]: 

1   Cp  TS  39 
Qz    2   6      3p   nu*  (        ' 

where  n  is  a  non-dimensional  parameter.   This  relationship  was 
used  in  both  rising  and  subsiding  air  representing  a  Fourier  expansion 
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of  the  heating  function.   The  factor  1/2  is  therefore  needed 
to  maintain  the  correct  amplitude  of  the  heating  perturbation. 
When  (2.6)  is  appled  at  level  2  the  result  is 

£■  (Y,  -  ¥,)  +  IKXVC4'1  ^3)-  V(Y1  -  V3)  -  Apo2a)2  =  0.   (2.8) 
Define  the  following  quantities 

^M  =  — "2—  •  <2'9) 

^T  =  Yl  2  *'  '  (2-10) 

which  implies  that 

^3   =  ^M  -  *T  •  (2'12) 

Here  H'T  is  proportional  to  the  layer  thickness  and  is  therefore 
a  measure  of  the  mean  temperature  in  the  layer.   Using  these 
definitions,  add  (2.4)  and  (2.5)  and  divide  the  result  by  two, 
obtaining 


-^   V2VM  +   |KXV¥M  •  V(V2VM)  +   |KXVyT  •  V(V2^T) 


+  3„4^  -  fer2-  =  0  (2.13) 

o  8x      2Ap 


For  the  second  forecast  equation,  subtract  (2.5)  from  (2.4) 
and  eliminate  w2  by  using  (2.8)  obtaining 
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■^(V2  -  y2)  ¥  T  +   |KX  VVM  •  V(V2  -  y2)^T  +  |K  X  V^  •  V(V2YM) 

9^T    fowe 
Po   3x     2Ap    u»  (2.14) 

where 

y2  =r^  •  (2.14a) 


These  are  the  prediction  equations  for  the  model. 

Although  the  normal  mode  method  gives  the  phase  speeds,  growth 
rates,  and  wave  structures  for  all  modes,  the  initial  value 
approach  gives  the  phase  speed,  growth  rate,  and  wave  structure 
of  only  the  most  unstable  mode.   This  approach  was  chosen 
because  only  information  about  the  most  unstable  mode  was 
needed. 

It  is  convenient  to  write  iK,  and  \pT    in  the  following  fo: 


nn 


\pM   =  E(y)  +  A(y,t)  cos  kx  +  B(y,t)  sin  kx ,  (2.15) 

tyj   -    F(y)  +  C(y,t)  cos  kx  +  D(y,t)  sin  kx.     .         (2.16) 

where  k  is  the  x  wave  number.   The  coefficients  A  through  D 
are  the  Fourier  amplitudes  of  the  disturbance  and  E  and  F  are 
the  mean  zonal  fields. 

These  expressions  for  \pM   and  \pT   are  substituted  into 
(2.13).   The  various  sine  and  cosine  terms  are  then  separated, 
and  all  products  of  the  quantities  A  through  D  neglected. 
Equating  the  coefficients  of  the  cosine  terms  gives 
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aT  ^ay7"  "  AJCJ   "  KL9y  ay^.  3y'8yT 


93ep      a_32 

9y3  D        8y3  u 


k2(f  B+f  D)]- 


3     Bk-    K(|-5- 
o  vdy2 


k2)(A-  C). 
(2.17) 


For  the  sine  terms  the  result  is 


,32B 


9E  92A    9F  92C    93E  A  ,  93F 


9t  ^9y2    DK    J      K  L"  3y  9y2    9y  dy2         9y3 


A  + 


9y3 


+  k2(|^A  +  14  C)]  +  3  A  k  -  K(|^5- 
v9y     3y  ;J     o       ^9y2 


k2) (B  -  D) 
(2.18) 


Repeat  the  procedure  for  equation  (2.14)  for  the  cosine  terms 
which  gives 


_9_  f92C 
9t  l9y2 


Ck2  -  CU2)  -  k[|£  |!?  +  |F  3fB  „  .|B 

'  L9y  9y2    9y  9y2    9y 


|-4  -  |f  D  (k2+  y2) 


9F 


9yB(k2"y2)    9y3-   9y3 


3'BD-  |i*B] 


3QDk  +  (K-  S)(|y2  -  k2)(A-  C). 


(2.19) 


The  results  for  the  sine  terms  are 


3      r92D 


9E 


9F 


ft    C|df  -Dk2    -   Dy2)    =   k  [|fC(k2    +    y2)     +    "  A  (k2    -     y2) 


M    9^C  9_F    9_^A         93F  93E  p, 

9y    9y2     "       9y    9y2         9y3  A        3y3  LJ 


where 


+    3QCk   +    (K    -    S)(|--   -k2)(B-  D), 


2RTC    LfJ 


%  v2 


(2.20) 


(2.20a) 
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and 

s  =  -f-  -ft- -g- -&<-£> '*•  (2-21> 

The  computational  equation  for  vertical  velocity  is  given 
below: 

9       dipT  di>T  di>r        vn 

M      =    _±_    [-JL   +    u     -3-i-   +    v     -^—   +    ^-]  ,  (2.22) 

2         Apo     L8t  m    3x  m    3y  2£0J '  v  J 


where 


um   =    "    iV  ~WL   '  (2'23) 

1         ^M 
vm   -^"-W-    '  (2'24) 
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III.   BOUNDARY  CONDITIONS  AND  FINITE 
DIFFERENCING  SCHEME 

The  finite  difference  scheme  is  illustrated  below  with 
a  sample  variable  M: 

f=  27T  CMi+l  -  Mi-l)  >  V'V 


l-T-   =  ur    (M-Al  -  2M.  +  M.  J  ,  (3.2) 

3y2     H2   *"  l+l        l      1-1/   '  v     J 


t^T  =    WT  t(M-o.9  "  2M-^i  +  M-)  "  (M-  -  2M.  ,  +  M.  9)] 
3y3    2H3  L v  i+2      i+l     i^    v  i      i-l     i-2^J 

(3.3) 
where  i  is  the  grid  index  and  H  is  the  distance  between  grid 
points . 

Centered  time  differences  were  used  for  all  quantities 
except  those  involving  friction  and  heating.   The  friction 
and  heating  terms  were  computed  at  time  (t  -  At) .   The  first 
step  in  all  cases  is  a  forward  time  step.   The  second  order 
equations  for  time  tendencies  were  solved  by  the  exact  method 
of  Richtmyer  [1957,  p.  101]. 

Yanai  and  Nitta  [1968a]  studied  the  problem  of  using 
finite  difference  equations  to  solve  dynamic  instability 
problems  of  non-divergent  barotropic  currents  with  boundaries. 
In  this  study  rigid  boundaries  are  placed  at  y  =  ±  W/2  where 
W  is  the  total  width  of  the  computational  field.   The  boundary 
conditions  were  A=B=C=D=Oaty=±  W/2. 
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Fig.  2.   Mean  wind  profile  U.  =  a.U   sech2  (y/L) 
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IV.   WIND  PROFILE  AND  INITIAL  CONDITIONS 


This  study  uses  the  following  zonal  wind  profile: 

U.  =  a.U   sech2  (y/L) ,  (4.1) 

j  j    o 

where  a.  can  be  ±  1  or  0.   The  barotropic  stability  of  this 
profile  has  been  analyzed  by  Lipps  [1962]  .   Figure  2  contains 
the  profile  (4.1)  for  a.  =  1. 

The  initial  conditions  are  given  below: 


A  =  sin  ny/W  (4. 2) 

B  =  0  (4.3) 

C  =  0  (4.4) 

D  =  0  (4.5) 

E  =  a,  +  a3 


F  =  a,  -  a3 


2 3.   foUQL  tanh  (y/L)  (4.6) 

j-^-   fQUoL  tanh  (y/L)  (4.7) 


Note  that  the  E  and  F  fields  are  independent  of  time. 
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V.   NON-DIMENSIONALIZATION  PROCEDURE 

In  order  to  simplify  the  interpretation  o£  the  results, 
the  equations  are  written  in  non-dimensional  form  as  follows 
x  =  Lx',  y  =  Ly%  and  t  =  L/t'.   The  stream  function  is 
written 


and 


^M  =  "  LUq/um  (y-)dy^  +  iPM(x',y  '  ,t ') 


yT  =  -  LUQ/uT  (y')dy*  +  *T(x*,y',t') 


(5.1) 


(5.2) 


when  these  expressions  are  substituted  into  (2.13)  and 
(2.14),  and  when  the  products  of  the  ^ ' s  are  dropped,  the 
non-dimensional  equations  can  be  written: 


-9 


d2u,  H 


<Tt   +  u4-^2^  +  Ce  -  3=B  j^1  +  rv2(^M  -  ♦  )  -  o  (5.3) 


&  +  u^Hv2  -  s)^T  ♦  (g  -  0)   ||I  +  (6-  r)v2C*M  -  ^ 

(5.4) 


=  0 


where 


-  i 


e  =  l/yz,  6  =  3  L ZU   ,  r  = 


and 


,  89   .   - 

<S    nmax  L  p0  g  u   3p  sin  2  a   Ae 


fpL  sin  2q   Ae  ,  '2 
U   2Ap/p4g  l2f0J 

v2 


Uo4  fo  T, 


2fo 


(5.5) 


(5.6) 


The  quantity  e  is  a  rotational  Froude  number,  6  is  a  non- 
dimensional  form  of  6,  T    is  the  ratio  of  the  square  root  of 
the  Ekman  number  [Greenspan,  1968]  to  the  Rossby  number,  and 


20 


6  is  the  heating  parameter.  The  primes  on  the  independent 
variables  have  been  dropped.   These  equations  are  equivalent 

to  the  set  2.17,  2.18,  2.19,  and  2.20  which  were  solved 
numerically. 
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VI.   COMPUTATIONAL  PROCEDURE 

The  following  constants  were  used  in  the  computations: 

p    =  100  centibars, 
*o 

L    =  400  kilometers, 

W    =  4000  kilometers, 

U    =10  meters  per  second, 
o  r 

f    =  varied  with  latitude, 
o 

a     =  3.0  meters  per  second  squared  per 
centibar  squared, 

3  =  varied  with  latitude, 

Ap  =  50  centibars, 

Ae  =  10  meters  per  second  squared, 

a  =22.5  degrees , 

p2  =  varied  with  latitude 

H  =50  kilometers, 

At  =1.0  hours , 

n     =  4,  6,  8 

max     '   ' 

8e73p=  -0.6  degrees  per  centibar. 
A  forecast  period  of  100  days  was  used  so  that  the  most  unstable 
mode  would  have  reached  the  maximum  growth  rate.   For  conven- 
ience, growth  rates  were  computed  over  the  last  two  days  of 
the  forecast  period.   This  was  done  by  assuming  the  amplitudes 
could  be  written  as 

x,  =  a0  exp  (n  UQ/L)  tx ,  (6.1) 

x2  =  aQ  exp  (n  UQ/L)  t2 ,  (6.2) 
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where  n  is  the  non-dimensional  growth  rate  and  a   is  a  constant 
By  forming  the  ratio  x2/x1,  the  growth  rate  can  be  shown  to  be 


t,  -  L  In  (x2/xi) 
n  '  U0  (t2  -  tj 


(6.3) 
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VII.   RESULTS 

A  series  of  experiments  were  conducted  over  a  wide  variety 
of  conditions.   The  more  important  experiments  are  given  in 
Table  I  as  a  function  of  ax  ,  a3  ,  3,  e,    6,  r,  n,  and  k%  where 
n  represents  the  maximum  growth  rate  for  each  set  of  conditions 
and  k'  is  the  corresponding  non-dimensional  wave  number  kL. 
When  some  of  these  parameters  were  computed  over  a  range  of 
values,  the  range  is  given. 

Experiment  1  was  conducted  to  determine  the  difference  in 
growth  rate  characteristics  between  easterly  and  westerly  jets. 
In  this  experiment  no  heating  or  friction  was  included  and  the 
mean  flow  was  the  same  at  both  levels.   In  order  to  obtain  a 
continuous  variation  between  the  easterly  and  westerly  flows, 
the  non-dimensional  parameter  3  (see  equation  5.5)  was  defined 
to  assume  the  sign  of  U.   Thus  $  is  positive  for  westerly  flow 
and  negative  for  easterly  flow.   In  the  remaining  experiments 
the  jet  will  be  easterly,  but  3  will  be  regarded  as  positive. 
For  each  3  the  maximum  growth  rate  and  the  corresponding  wave 
number  were  determined.   Figure  3  contains  these  quantities 
as  a  function  of  3;  the  maximum  growth  rate  of  n  =  0.1888  and 
occurs  at  3  =  -0.352.   The  corresponding  wave  number  is  k'  = 
0.99.   This  value  of  k'  is  close  to  the  minimum  k"  which  occurs 
when  3  =  -0.112.   Fig.  3  shows  clearly  that  the  easterly  jet 
is  more  unstable  than  the  westerly  jet.   In  particular,  the 
growth  rates  for  the  easterlies  over  a  range  of  beta  are 
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greater  than  the  value  for  beta  =  0.   This  is  in  contrast  to 
the  normal  notion  that  the  presence  of  beta  has  a  stabilizing 
effect  for  instabilities. 

Experiment  2  investigated  the  downward  propagation  of 
energy  from  a  barotropically  unstable  disturbance  as  a  func- 
tion of  the  latitudinal  location  of  the  jet.   An  easterly  jet 
was  introduced  at  the  upper  level  and  the  lower  level  mean 
flow  was  set  to  zero.   Surface  friction  was  included,  but 
heating  was  neglected.   A  measure  of  the  downward  propagation 

is  the  ratio  R  =  I  ¥,  I    /l^il     where  these  amplitudes  are 

1  3 'max  '  i 'max  r 

obtained  from  the  exponentially  unstable  solution.   Figure  4 

contains  isolines  of  R  as  a  function  of  k'  and  latitude.   The 

curve  n    shows  the  wave  number  which  has  the  maximum  growth 
max  b 

rate  at  each  latitude.   The  R  value  decreases  as  one  moves 
from  high  to  low  latitude  as  predicted  by  Charney  [1963]  . 
The  R  value  also  decreases  with  increasing  wave  number  at 
each  latitude.   This  figure  shows  the  dependence  on  epsilon. 
Experiments  3,  4,  and  5  all  contain  heating  and  were  con- 
ducted with  a  fixed  latitude  of  10  degrees,  otherwise  they 
are  the  same  as  experiment  2.   In  these  experiments  a  weighting 
factor  was  applied  to  the  heating  following  Bates  [1970]  and 
Yamasaki  and  Wada  [1972] .   The  following  expression  was  used 

n  =  n    /cosh2y/L.  (7.1) 

•    'max      } '  K        J 

This  form  implies  that  moisture  is  available  for  the  heating 
in  the  region  |y|<  L.  If  this  was  not  done,  CISK  heating  would 
occur  in  regions  away  from  the  basic  jet.   In  experiment  3 
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where  n     =  4  the  non-dimensional  heating  parameter  is 
lmax  to  r 

6  =  0.165.   Figure  5  shows  the  growth  rate  (n)  and  the  ampli- 
tude ratio  (R)  as  a  function  of  the  non-dimensional  wave  number 
k" .   Note  that  the  maximum  growth  rate  of  n  =  0.187  occurs  at 
k'  =  1.0.   These  two  values  are  very  close  to  the  values  ob- 
tained without  heating.   The  R  is  fairly  large  at  the  lower 
values  of  k'  with  a  maximum  at  0.5.   R  becomes  very  small  in 
the  area  of  maximum  growth.   Apparently  the  barotropic  insta- 
bility effects  dominate  in  the  areas  of  maximum  growth  while 
the  heating  is  relatively  more  important  for  k'  £  0.5  as 
evidenced  by  the  larger  values  of  R.   This  effect  is  even 
more  pronounced  in  experiment  4. 

In  experiment  4  where  n     =6  the  non-dimensional  heating 
^  'max  & 

parameter  is  6  =  0.247.   Figure  6  shows  that  the  maximum  growth 

rate  of  n  =  0.186  occurs  at  k^  =  1.0  just  as  in  experiment  3. 

As  k'  varies  away  from  1.0  to  either  extreme  the  growth  rate 

approaches  0.167.   The  R  is  larger  in  both  of  these  regions 

where  the  largest  value  occurs  for  the  largest  value  of  k'. 

A  secondary  maximum  in  R  occurs  at  k'  =  0.6.   This  indicates 

that  heating  dominates  the  solution  in  the  regions  which  have 

the  larger  values  of  R. 

In  experiment  5  where  n    =8  the  non-dimensional  heating 
1  'max  & 

parameter  is  6  =  0.330.   Figure  7  shows  that  the  growth  rate 
is  fairly  constant  over  the  k'  considered.   This  is  similar 
to  the  results  obtained  by  Williams  and  Robertson  [1973]  who 
acquired  essentially  a  constant  growth  rate  with  a  hyperbolic 
tangent  wind  profile  using  larger  values  of  6. 
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GROWTH  RATE  n 


Fig.  j> .   The  growth  rate  n  (solid  line)  and  the  corre- 
sponding wave  number  k/ (dashed  line)  as  functions 
of  $  for  experiment  1.   Westerly  flow  occurs  for 
3  >  0  and  easterly  flow  for  3  <  0. 
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EPSILON 


+  0.6 


Fig.  4.   Isolines  of  R  =  |y,|    / I V , |     (curved 

b  '  3  ■ max  '  i ' max 

lines)  as  functions  of  k",  latitude,  and 
eps  ilon. 
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Fig.  5. 


Growth  rate  (dashed  line)  and  amplitude  ratio 
(solid  line)  as  a  function  of  wave  number  for 
experiment  3. 


30 


GROWTH  RATE  n 


AMPLITUDE  RATIO  R 

Fig.  6.  Growth  rate  (dashed  line)  and  amplitude  ratio 
(solid  line)  as  a  function  of  wave  number  for 
experiment  4. 
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Fig.  7.  Growth  rate  (dashed  line)  and  amplitude  ratio 
(solid  line)  as  a  function  of  wave  number  for 
experiment  5. 
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VIII.   CONCLUSIONS 

Recent  studies  suggest  that  disturbances  may  form  through 
barotropic  instability  in  the  easterly  jet  south  of  the 
Tibetan  High  at  200-mb.   In  this  thesis  the  instability  of  an 
analytical  jet  profile  was  investigated  with  a  two-level 
quasi-geostrophic  model.   The  downward  propagation  of  energy 
from  a  barotropically  unstable  upper  level  disturbance  was 
investigated  to  see  if  it  could  initiate  lower  tropospheric 
disturbances . 

The  hyperbolic  secant  squared  wind  profile  was  used  in 
all  of  the  numerical  experiments.   The  first  experiment,  which 
contained  no  heating  or  friction,  compared  the  growth  rates 
with  easterly  and  westerly  flow.   The  easterly  jet  was  the 
more  unstable  and  its  maximum  growth  rate  was  well  above  the 
growth  rate  for  3=0.   The  energy  propagation  was  studied 
by  introducing  the  jet  at  the  upper  level  with  no  mean  flow 
below.   The  ratio  (R)  of  the  disturbance  amplitude  at  the 
lower  level  to  the  disturbance  amplitude  at  the  upper  level 
was  computed.   With  surface  friction  included,  R  decreased 
from  high  latitudes  to  low  latitudes;  at  each  latitude  R 
decreased  with  increasing  wave  number.  When  CISK  type  heating 
was  introduced  with  a  large  enough  heating  coefficient,  the 
growth  rate  became  independent  of  wave  number.  In  this  case 
R  had  a  finite  value  and  it  increased  with  increasing  wave 
number.   For  smaller  values  of  the  heating  coefficient  the 
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barotropic  instability  controls  the  growth  rate  and  R  remains 
very  small. 

The  result  that  the  easterly  jet  is  more  unstable  than 
the  westerly  jet  supports  the  view  that  barotropic  instabil- 
ity occurs  south  of  the  Tibetan  High.   The  small  value  of  the 
amplitude  ratio  at  low  latitudes  agrees  with  Charney's  [1963] 
conclusion  that  the  vertical  energy  propagation  is  small. 
The  inclusion  of  CISK  heating  did  increase  the  low  level 
amplitude,  but  this  could  not  occur  unless  a  finite  amplitude 
were  present  at  the  lower  level.   Since  these  results  were 
obtained  with  a  quasi-geostrophic  model  with  only  two  levels, 
it  is  important  to  carry  out  further  studies  using  more 
levels  and  also  allowing  for  gravity  waves.   This  would 
require  a  multi-level  primitive  equation  model  in  either 
spherical  coordinates  or  on  an  equatorial  6-plane.   A  more 
sophisticated  heating  parameterization,  such  as  that  by 
Arakawa  and  Schubert  [1974]  ,  should  be  used. 
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